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Abstract 

We present analytical expressions and numerical illustrations for the ground-state 
density distribution of an ideal gas of spin-polarized fermions moving in two dimen- 
sions and driven to rotate in a harmonic well of circular or elliptical shape. We show 
that with suitable choices of the strength of the Lorentz force for charged fermions, 
or of the rotational frequency for neutral fermions, the density of states can be 
tuned as a function of the angular momentum so as to display a prominent shell 
structure in the spatial density profile of the gas. We also show how this feature of 
the density profile is revealed in the static structure factor determining the elastic 
light scattering spectrum of the gas. 
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1 Introduction 



Ultracold Fermi gases of alkali atoms such as 40 K or 6 Li in harmonic traps 
are quantum systems that are experimentally accessible by the techniques of 
atom trapping and cooling [1,2]. The high purity and the low temperature of 
the samples and the high resolution of the detection techniques make these 
systems ideal candidates for the study on the mesoscopic scale of single-level 
quantum properties such as shell structures in the particle density profiles 
[3-7]. In the experiments the trapped atomic gas can be fully spin-polarized 
and the strength and the anisotropy of the trap can be tuned to reach quasi- 
onedimensional (ID) or quasi-t wo dimensional (2D) configurations. Quantum 
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effects in the equilibrium profiles can be greatly enhanced by varying the 
anisotropy of the confinement [6]. Here we show how to enhance quantum 
effects on 2D equilibrium density profiles by putting the gas into rotation. 

In a wholly different physical realm electrons in a quantum dot realize a Fermi 
gas moving in a plane under the effect of a harmonic potential. A quantum 
dot is often described as an artificial atom, in which a confinement potential 
replaces the attractive potential of the nucleus and the increased length scale 
allows access to experiments that are not practicable on real atoms [8-10]. Spin 
effects are also accessible and complete spin polarization is easily achieved in 
a small dot. Confinement by a circularly symmetric harmonic potential is of 
special interest, since in this case the energy needed to add electrons to the dot 
reveals a shell structure [11]. The effects of a magnetic field on the electronic 
properties of circular quantum dots have been studied by numerical diagonal- 
ization and by spin-density functional methods (see Ref. [12] and references 
therein). These approaches have also been applied to rectangular dot struc- 
tures, when the effective lateral confinement is described as elliptical. In such 
a confinement the single-particle energy levels as well as the energy spectra 
for two-electron systems are exactly known as functions of the magnetic field 
[13,14]. 

A Fermi gas of spin-polarized charged particles in a uniform magnetic field, 
under conditions such that the Coulomb interactions can be neglected as speci- 
fied below, can be mapped into a rotating Fermi gas of neutral atomic particles 
in a state of complete spin polarization, where the atom-atom interactions are 
negligible on account of the Pauli principle suppressing s-wave scattering. In 
this work we examine the particle density profiles of the ideal system that 
we have just introduced. For definiteness we shall refer throughout the paper 
to the ideal 2D Fermi gas of charged particles in a uniform magnetic field, 
although some parts of the discussion (and in particular the analysis of the 
elastic light scattering spectrum given in Sec. 4) are specifically aimed at an 
atomic Fermi gas in uniform rotation. 

An ideal 2D Fermi gas of charged particles in a uniform magnetic field has 
recently been studied by van Zyl and Hutchinson [15], who have used an in- 
verse Laplace transform method and an expansion in Laguerre polynomials 
to calculate its thermodynamic properties. Their method can be implemented 
exactly in the case of a uniform gas, whereas recourse to a local-density ap- 
proximation is needed for an inhomogeneous gas on account of the broken 
symmetry induced by the concurrent effects of the harmonic confinement and 
of the magnetic field. However, the single-particle Hamiltonian of the ideal 
gas can be diagonalized on the basis of left and right circular or elliptical 
quanta, bringing it to the form of the Hamiltonian of two independent har- 
monic oscillators. The eigenfunctions can thus be written in terms of Hermite 
polynomials, as we do in this work. Practical applications of our approach are 
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limited to ultracold Fermi gases consisting of a restricted number of particles, 
so that the number of Hermite polynomials that need explicit numerical calcu- 
lation remains limited. We shall give below numerical illustrations for systems 
of 10 and 106 fermions at zero temperature. 



The paper is organized as follows. In Sec. 2 we set out the single-particle 
Hamiltonian of the system in a generic gauge which allows its diagonalization 
on the basis of left and right-handed circular or elliptical quanta. In Sec. 3 
we express the spatial density profile of the gas in terms of Hermite polyno- 
mials and display several configurations for various values of the strength of 
the magnetic field and of the anisotropy parameters of the confinement. We 
show that an impressive shell structure arises in the profile whenever there 
is an excess of density of states with zero angular momentum. In Sec. 4 we 
calculate the elastic contribution to the static structure factor of the gas, in 
order to show how the shell structure is reflected in the elastic light scattering 
spectrum. Finally Sec. 5 offers some concluding remarks. 



2 The model 



We consider a 2D spin-polarized Fermi gas at zero temperature, made up of 
iV non- interacting fermions of charge e and mass uif subject to a uniform 
magnetic field B = V x A perpendicular to the plane of confinement (the 
{x,y} plane, say). The gas is immersed in a uniformly charged background 
medium so that the total system is neutral. Screening effects cut the range of 
the Coulomb interactions and to a first approximation these can be neglected 
on account of the Pauli principle keeping fermions with parallel spins far apart 
from each other. More generally, the Coulomb interactions between electrons 
in a quantum dot are negligible when the harmonic oscillator length charac- 
terizing the confinement is small compared with the effective Bohr radius. 

As already noted in Sec. I, this system can be mapped into a spin-polarized 
gas of N neutral fermions confined in a rotating trap. The map is guaranteed 
by a correspondence between the Maxwell equations and the Navier-Stokes 
equations in the absence of viscosity, as demonstrated by Marmanis [16]. In the 
rotating system the velocity field v of the fermions and the vorticity w = Vxv 
play the role of the vector potential A and of the magnetic field B, respectively. 
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2.1 Circular trap 



It is useful to present first the case of an isotropic harmonic potential. The 
single-particle Hamiltonian for non-interacting fermions is 

In Eq. (1) ujq is the trap frequency, uj l = eB/(2cm,f) is the Larmor fre- 
quency, and we have chosen the symmetric gauge for the vector potential 
A = (— By/2, Bx/2, 0), as is suited to the circular symmetry of the confine- 
ment [17]. 

The Hamiltonian (1) is separable [18] upon introducing first the creation and 
destruction operators a x y and a x>y for a quantum of oscillation in the x or 
y direction, defined according to a x = [(^/m/fi)" 1 / 2 ^ + i{hmfVL)~ l l 2 p x ]/ \fz~ 
where VL = [u 2 + u 2 L ) l l 2 , and then the operators a d = (a x — ia y )/\^2 and 
a g = (o-x + ia y )/\/2. The transformed Hamiltonian is 

H = hw d (a d a d + + hw g (aj,a 9 + (2) 

where u d = Q + ul and ui g = Q — ojl- 

The operator a d (a g ) can be interpreted as the destruction operator of a right 
(left) circular quantum, the number of such quanta being given by h d = a\a d 
{fig = a g a g ). The quantity h(n d — h g ) gives the eigenvalues of the angular 
momentum operator L z = xp y — yp x . 



2.2 Elliptical trap 



In the case of anisotropic 2D confinement neither the symmetric gauge nor 
the Landau gauge are useful. We write the single-particle Hamiltonian in the 
general form 

H = ^{p x -aeBy/c) 2 + ^{p y + {l-a)eBx/c) 2 + lm f uj 2 {x 2 + X 2 y 2 ) (3) 
Zm f Zm f Z 

with < a < 1 and A the anisotropy parameter. Equation (3) can be rewritten 

as 

1 1 

H = 7^(PI + P 2 y) + 2 m f^l^ 2 + ^If) + 2uj l [(1 - a)xp v - ap x y] (4) 
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where we have set oo 2 = u 2 + 4(l — a) 2 u)1 and u> 2 = X 2 uj 2 + 4a 2 uj1. By imposing 
the condition (1 — a)/a = u x (a) / u y (a) , which has the solution 

a = A/(l + A), (5) 

and using a x = [{h/mfUJ x )~ l l 2 x + i(hmfU x )~ 1 ^ 2 p x ]/y/2 etcetera, Eq. (4) be- 
comes 



H = ^h (u x + toy) (d x d x + ciydy + 1) + (u x - u y ) (dld x - a^a^) 

+ 2 ^ ihw L \a x al - a x a y \ . (6) 
1 + A L J 

For circular confinement (A = 1) the a-gauge given by the relation in Eq. 
(5) reduces to the symmetric gauge, while in the limit of one-dimensional 
confinement (A = or oo) it reduces to the Landau gauge. More generally, 
in the mapping with a rotating system the a-gauge corresponds to a velocity 
field which follows the symmetry of the confinement. 

The Hamiltonian in Eq. (6) is separable in terms of the operators 

dd = cos 9 a x — i sin 6 a y 
ci g = sin 9 d x + i cos 9 d y 

if the angle 9 satisfies the relation tan(26 ) ) = 2ujl/{uj x — u y ), with lj l = 
2y/Xu L /(l + X). In the new basis the Hamiltonian can be written as in Eq. (2) 
with 

UJ d ,g = ^K + U y ) ± + j(uJx-Wy) 2 (8) 

The operators and d g are the destruction operators for a right-handed 
and a left-handed elliptical quantum, but (if 9 ^ 7r/4) the quantum number 
m = n d — n g is no longer proportional to the angular momentum L z . 



3 Spatial density profiles 



The eigenstate \xn d ,n g ) corresponding to the occupation numbers {rid,n g } can 
be obtained by recursively applying the operators a d and d g to the ground 
state |xo,o)- Using Eq. (7) we can write 

\Xn d ,n g ) = —j=L=(cos9al + i sin 6»aJ)™ d (sin 6»4 - i cos6»aj) ns |xo,o) , (9) 
^n d \n g \ 
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where 



(x,v\Xofi) = n~ 1/2 ex P (-x 2 /2l 2 x )eM-y 2 /2ll) (10) 

with l XtV = ^Jh/rrifUJ Xt y. Series expansion of Eq. (9) and use of the expression 
of the eigenstates of the harmonic oscillator in terms of Hermite polynomials 
lead to the result 



Xn d ,n g { x }y) — i t ; ^ — 2-~i • 



n d l n g ! 



— (-1)™ 9 ~ 

(i)na + n g -k- j{sin9) n d - j+ k {cos9) n g -k + j H ^ 

(11) 

Finally, the density profile n(x,y) can be expressed in the form 



n 



N s N d 

fay) = J2 \Xn d ,n g (x,y)\ 2 . 

n g =0 n d =0 



Here, 



N d = Int 



1\ LU n 



— I n J — 

hu d V 9 1) u d 2. 
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12) 



(13) 



is the highest allowed value for the quantum number n d at given values of the 
Fermi energy E F and of the quantum number n g , and 



N g = Int 



u d 



fojJg 2lJg 



(14) 



is the highest allowed value for the quantum number n g at given Ep. Of course, 
the Fermi energy is determined by the number of fermions in the trap. 

Figure 1 shows the particle density profile for 10 spin-polarized fermions in an 
isotropic trap (A = 1) at various strengths of the applied magnetic field, and 
Fig. 2 reports the corresponding density of occupied states (dos). In zero field 
(Figs. 1(a) and 2(a)) the total angular momentum is zero and the left and 
right states are equally populated. At a field strength such that the Larmor 
frequency is equal to the trap frequency (Figs. 1(b) and 2(b)) there still is 
an excess population of the state m = and a central peak in the density 
profile. At still higher fields (Figs. 1(c) and 2(c)) the dos becomes constant 
and the occupation of two states with opposite angular momentum (m = ±1) 
produces a central hole in the density profile. Finally, at very large fields (Figs. 
1(d) and 2(d)) the Landau quantization limit is being recovered, with a filling 
factor equal to unity in a dos which is constant for m < 0. This is spatially the 
most compact distribution of the droplet of spin-polarized fermions, which is 
commonly denoted as MDD [19]. 

Similar density distributions can be obtained for higher numbers of fermions. 
In Fig. 3 we report as an example the spatial distribution of 106 fermions at 
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uol = 5u;o, showing as in Fig. 1(b) a central peak from an excess dos at m — 0. 
On a slight increase in the field strength, to ujl > 5.1uj , the filling factor 
becomes unity and the density profile takes a MDD configuration similar to 
that shown in Fig. 1(d). 

In an anisotropic trap the overall elliptical shape of the fermion cloud can still 
sustain a shell structure, which becomes progressively smoothen as the field is 
increased. In Fig. 4 we show the case of 10 fermions in a trap with anisotropy 
parameter A = 4. For uj l = uj (Fig. 4(b)) the dos is similar to that shown for 
a circular trap in Fig. 2(b), but the value of m is not related to the angular 
momentum in the anisotropic trap and a central peak in the density profile 
is consequently missing. Instead, raising the field produces a softening of the 
oscillations of the profile along the x axis and some enhancement of those 
along the y axis. This effect is more pronounced at wl = 1.8u;o (Fig. 4(c)), 
where a central hole in the profile is accompanied by two bumps along the y 
axis. For this value of ujl the dos is similar to that shown for a circular trap 
in Fig. 2(c). Again, all shell structures disappear and an MDD-like profile is 
obtained when the filling factor becomes unity (Fig. 4(d)). 

Finally, Fig. 5 reports the density profile for 10 fermions in a strongly anisotropic 
trap (A = 10). In zero field (Fig. 5(a)) the system is quasi-onedimensional: 
only the ground state is occupied in the y direction and density oscillations 
are present only along the x axis. These become progressively smoothen as 
the field is applied (Fig. 5(b)) and increased (Fig. 5(c)). 



4 Angular distribution of scattered light 

The theory of light scattering from a confined cloud of spin-polarized fermionic 
atoms has been presented in Ref. [20]. Here we discuss an experiment of elastic 
scattering in which the incident beam is propagating along the z direction, 
orthogonally to the {x, y} plane of confinement of the gas. 

The positive-frequency component of the incident electric field is 

E^(z,t) = -Eiexp[i(k L z-ujt)] (15) 

with E and e being the field amplitude and polarization vector, and k^ and 
uj its wave number and frequency. The angular distribution of the scattered 
light can be decomposed into two components: an elastic term originating from 
diffraction by a finite object of given optical density profile, and an inelastic 
term determined by excitations in the quantum fluid. Only the elastic term 
remains if the light frequency is far-off-resonance from the internal states of the 
fluid. This term is determined by the elastic static structure factor S e (k x , k y ), 
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which is defined by 



S e (k x ,k y ) = J dx J dye ikxX e ikyV n(x,y) 
where k = (k x , k y ) is the wave vector transfer in the scattering process. 



(16) 



We have evaluated the elastic contribution to the static structure factor for 
the density distribution of 106 fermions shown in Fig. 3 and for a MDD con- 
figuration found with the same number of fermions at Ul = 5.1uj . These two 
configurations have essentially the same size and thus show almost identical 
central peaks in their diffraction pattern (see the inset in Fig. 6). According 
to Babinet's principle this peak is given by the optical diffraction pattern of 
a circular aperture. 

The central part of the density profile shown in Fig. 3 consists of a peak of 
diameter ~ 2l x surrounded by a ring of diameter ~ 4l x . As noted in Sec. 
3, this double structure arises from an excess of occupied states with zero 
angular momentum and its signal appears in the diffraction pattern in the 
range 1.5 < k x l x < 3. In the same range the scattered intensity from the flat 
MDD configuration is much lower, as is seen in the main frame of Fig. 6. 



5 Summary and conclusions 



In summary, we have evaluated the particle density profiles of an ideal gas 
of spin-polarized fermions confined in 2D inside a circular or elliptical har- 
monic trap as functions of an applied magnetic field or of the trap rotational 
frequency. We have introduced a generic a-gauge in which the eigenstates of 
the Hamiltonian can be expressed in terms of right-handed and left-handed 
elliptical quanta and written in terms of Hermite polynomials. The a-gauge 
reduces to the symmetric gauge in a circular trap and to the Landau gauge in 
a strongly elongated trap. 

In the case of a circular trap we have correlated the shell structures in the den- 
sity profile of the gas with the density of occupied angular-momentum states 
and suggested that these structures may become observable in an elastic light 
scattering experiment from an atomic Fermi gas inside a rotating pancake- 
shaped trap. Other peculiar shell structures have been illustrated in the case 
of an elliptical trap. It remains a challenge to provide an exact treatment for 
a gas consisting of very large numbers of particles at finite temperature. 
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Fig. 1. Density profile for 10 fermions in a circular harmonic trap at (a) wx, = 0, (b) 
ujl = wo, (c) ujl = 1.18cjo, and (d) lol = Wojq. 



(a) 



(b) 



(c) 



(d) 
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Fig. 2. Density of occupied states (dos) as a function of the quantum number M 
rid ~ n g f° r each of the configurations shown in Fig. 1. 




Fig. 3. Density profile (left) and dos (right) for 106 fermions in a circular harmonic 
trap at lol = 5wq- 
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Fig. 4. Density profile for 10 fermions in an elliptical harmonic trap with anisotropy 
parameter A = 4 at (a) u>l = 0, (b) u>l = u>q, (c) ujl = 1.8wq, and (d) ojl = 10u>o- 




Fig. 5. Density profile for 10 fermions in an elliptical harmonic trap with strong 
anisotropy (A = 10) at (a) u>l = 0, (b) ujl = wq, and (c) u>l = Wuq. 




Fig. 6. Elastic contribution to the static structure factor S e (k x ,k y ) corresponding 
to the spatial density distribution shown in Fig. 3 (continuous line) and to an 
MDD distribution of the same number of fermions (dashed line), as a function of 
(k x l x + kyl y )/2 at k x l x = k y l y . The inset shows the central peak of the two patterns. 
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